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Abstract. The application of a classical approach to various quantum problems - the 
secular perturbation approach to quantization of a hydrogen atom in external fields 
and a helium atom, the adiabatic switching method for calculation of a semiclassical 
spectrum of hydrogen atom in crossed electric and magnetic helds, a spontaneous decay 
of excited states of a hydrogen atom, Gutzwiller's approach to Stark problem, long- 
lived excited states of a helium atom recently discovered with the help of Poincare 
section, inelastic transitions in slow and fast electron-atom and ion-atom collisions - is 
reviewed. Further, a classical representation in quantum theory is discussed. In this 
representation the quantum states are treating as an ensemble of classical states. This 
approach opens the way to an accurate description of the initial and final states in 
classical trajectory Monte Carlo (CTMC) method and a purely classical explanation 
of tunneling phenomenon. The general aspects of the structure of the semiclassical 
series such as renormgroup symmetry, criterion of accuracy and so on are reviewed 
as well. In conclusion, the relation between quantum theory, classical physics and 
measurement is discussed. 
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1. Introduction 



Classical physics is the foundation-stone in understanding of a microcosm, since 
all experimental devices are designed on classical principles. A classical approach 
can be used as an approximation in quantum physics also providing, sometimes, a 
more adequate description of dynamics than standard quantum approximations. The 
presented treatment involves all necessary quantum properties and the word 'classical' 
in the title of this paper is used just to emphasize that this review (beside Sec. 5) is based 
on the analysis of classical trajectories - not on the Schrodinger equation in semiclassical 
approximation. Besides, the classical approach to the real physical problems, for which 
an accurate asymptotic description can be developed, is presented, and original papers 
are cited only. The abstract problems, such as billiard models, Feigenbaum Universality 
and etc., are not discussed here since they are rather part of mathematics than physics. 

The classical approach plays a fundamental role in quantum mechanics. Our 
understanding of any result in quantum theory relies upon classical language. For 
instance, since the degree of freedom 'spin' has no analogue in classical mechanics, we 
have no idea how it looks like. The classical approach has also an advantage over 
different quantum approximations because of a more adequate description of dynamics 
of the system. 

The semiclassical quantization conditions are formulated in the configuration space: 
p(s)ds = ir(k + oti + a 2 )h , k = 0,1,2,... (1.1) 

_ i 

where Si and s 2 are caustics of a classical motion along the variable s of the coordinate 
system where the Schodinger equation can be separated, and a\ and a 2 are purely 
quantum phase shifts (so-called Morse indices) originating from the caustics; in the case 
of turning point a — |, in the case of Coulomb singularity a = — \ and in the case 
of the rotational motion a = 0. In contrast to classical mechanics, the quantization 
condition (jl.ip is not invariant with respect to canonical transformations pQ. For 
instance, in action-angle variables {/, if} obtained from the set {p q , q} by a canonical 
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transformation, the motion is always a uniform rotation (I = const, if = ut) and we 
lose the information about the type of caustics (turning points) which play a key role 
in quantization conditions. 

For the isotropic potential V(r) the problem is separable in the spherical coordinate 
{r, if} and Hamiltonian takes the form 



H 



2m e 



Pr 2 + - 2 {P* 2 + ^L)}+V(r), (1.2) 

r z \ sin ii ' J 



r 2 ^ sin $ ' 



where m e is the mass of the particle, p r = m e r, p$ = m e r 2 $ and p v = m e r 2 sin 2 $(p 
are canonically conjugate momenta to coordinates r, $ and if, respectively. Since 
particle rotates uniformly along ip, p v is conserved, caustics are absent and quantization 
condition (jl.ip gives[j]: 

p v = mh , m = 0, ±1, ±2, (1.3) 

The motion along $ variable is oscillation between two turning points i9i = arcsin(m^/L) 
and $2 = 7T — aicsm(mh / L) , where L is an angular momentum which is obtained from 
the quantization condition 



*2 L2 _ ™^}_ M = n + \/2)h , l=\m\,\m\ + l,..., (1.4) 
i V sin v 

as L = (I + 1/2)% beside one exception -1 = 0. In this case we should accept value 
L = which follows from quantum treatment^. At last, the energy level E Hr i for I ^ 
is determined from the radial quantization condition 




E nrl - V(r) - -^]dr = (n r + l/2)h , n r = 0, 1, 2, .... (1.5) 

The quantization condition (II. 5ft gives correct value of energies even for the low lying 
levels since in this case the actual potential can be approximated by a parabolic well 
for which semiclassical and quantum spectra coincide. In the case 1 = (L = 0), the 
quantization condition (II. 5p depends on the type of inner caustic r\. For Coulomb 
attraction a\ = — 1/4 and it compensates the contribution from the outer turning point 
«2 = 1/4 providing correct energy value. 

Another general aspect is the use of the asymptotic techniques. The semiclassical 
approximation or the perturbation theory is the asymptotic series over a small parameter 
^ : Hn a n^ n - In the first case A = h 2 , and in the second A is the amplitude of 
perturbation. Originally the asymptotic term {a n X n } is the estimation of the accuracy 
in the (n — l)-st order. At the beginning the terms a n \ n decrease with increasing n 
up to some index n = n m j n (A). After that the terms start to increase and we lose the 
accuracy. The term a n X n can be considered as a correction only if the next term is 
smaller. Since the larger A the smaller n min (X) is, in application, the leading term gives 

| In (jl.ip integration corresponds to the half-period, i.e. from ip = up to ip = it 

§ The semiclassical quantization rule is asymptote at I >> 1, and formally this case is beyond the 

validity of semiclassical approximation. 
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the correct result in the widest interval of A and we do not need any corrections which 
gradually lose the meaning of correction with increasing A. 

Atomic units (m e = e = h = 1) are used throughout the review, unless otherwise 
explicitly indicated. 

2. Energy spectrum and resonances 

2.1. Secular perturbation theory 

In the secular perturbation theory the Hamiltonian has the form H(p, r) = H (p, r) + 
XV (r), where XV (r) is perturbation. It is assumed that the unperturbed Hamiltonian 
Ho(p,r) is separable. In the iV dimensional case the unperturbed system has N 
independent integrals of motion Aj (i = 1,2, ...,N), including the energy. Under 
weak perturbation these integrals of motion begin to change slowly in time. Then, 
the equation of motion for Aj(t) can be averaged over an unperturbed motion (at fixed 
values of A;) from which corrected integrals of motion and quantization conditions are 
obtained. 

Usually, the secular perturbation theory is formulated in action-angle variables 
{/, (p} which are obtained from the set {p q , q} by a canonical transformation (see, 
e.g., [2]). However, these variables are not appropriate for quantization conditions, 
because after the canonical transformation the motion is always a uniform rotation 
and the information about the type of caustics and Morse indices is lost[JJ]. 

2.1.1. Hydrogen atom in crossed electric and magnetic fields. When two charged 
particles are moving in uniform electric F and magnetic B fields the problem of 
separation of the center of mass motion becomes nontrivial (see, e.g., [4]). The 
Hamiltonian of two particles with charges Z\ = e and Z 2 = —e in external electric 
and magnetic fields has the form 

H = 57> + l Al f + i( p2 - l M + eF(ri - r2) - 7 (2 ' 1} 

where r = r\ — r 2 , c is the light velocity, mj, r i; p^, and Aj are the mass, the 
radius- vector, the momentum and the vector-potential for i-th particle (i = 1,2). The 
momentum is connected with the velocity Vi by the relation pi = Vi — (e/c)Aj. The 
gauge for the vector-potential is chosen in the form Ai = [B x T*j]/2, where squared 
brackets [a x b] denote the vector multiplication of two vectors a and b. In this case, 
instead of the total momentum the vector P = pi + p 2 — (e/2c)[B x r] is conserved. 
After separation of 'the center of mass motion' the Hamiltonian for a relative motion of 

|| For instance, in [3] the action-angular variables are used to obtain the semiclassical energy spectrum 
in the quadratic Zeeman effect, but semiclassical quantization condition (3.8) in this paper has wrong 
Morse indices: a\ — 012 = instead of «i = aa = \ (see equation (12.19)) in Sec. 2.1.2) and leads to an 
incorrect energy spectrum (compare table 3 in [3] with table 1 in Sec. 2.1.2). 
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particles takes the form 

H = i p2 - i~J- B ■ t- x p » + 8^ |B x r| + « F + F <"> ' ') " 7< (2 - 2) 

where M = mi + m2 is the total mass, m = m\m2j{jn\ + m-i) is the reduced 
mass, \i = m\mil (mi — 777,2) and (a • b) denotes the scalar product of two vectors 
a and b. The separation constant P gives rise to an additional effective electric field 
F e ff = [P x B]/Mc. The appearance of the effective field F e ff is the trace of gauge 
invariance in the two body system; it reflects the uniform character of the homogeneous 
magnetic field. 

The specific feature of hydrogen-like initial Hamiltonian Hq is its huge degeneracy. 
In classical mechanics this degeneracy is manifested as a closed elliptic orbit of electron. 
In this case, we have only one quantization condition along this orbit, and only one 
quantum number is well defined. The missed quantization conditions are determined 
by the perturbation. For instance, for the spherical perturbation, the quantization 
conditions are written in the spherical coordinates, but if the perturbation is the uniform 
electric field, the quantization conditions are written in the parabolic coordinates. 

The semiclassical quantization of the hydrogen atom in weak electric and magnetic 
fields at arbitrary orientations of F and B was carried out by Epstein in 1923 [5J. In 
this case, the equation of motion for the electron is 

% = -- 3 -F+kBxv], (2.3) 
at r 6 c 

where v is the velocity of the electron; in the nonrelativistic case it is connected with 
momentum as p = m e v. The unperturbed Kepler elliptic orbit is specified by angular 
momentum L = [r x p] and Runge-Lenz vector A = [L x jj] + r/r, and can be presented 
in terms of L and A as 

r(t) = a[cos£(t) - e]j + ay/l - e 2 singft) ^ (2-4) 



where a is the semimajor axis of the ellipse, e = y 1 — L 2 /n 2 is the eccentricity and £ is 
the Kepler anomaly (or 'elliptic time'), which is connected with the actual time by the 
relation t = y/a(£ — esin£). Under influence of perturbation L and A change slowly in 
time and the equations of motion averaged over Kepler period read 

dL 3 Or . t 1 r - T t g?t4. 3 r _ -, X r . - , _ _ , 

— = - n 2 F x A + — B x L , — = - \F x L + — \B x A . 2.5 

dt 2 1 J 2c 1 J ' dt 2 [ J 2c 1 

Now let us introduce instead of L and A the new variables J12 = (L ± nA)/2 which 
are subject to the relation 

1 n 2 

J - = -m = t (2 - 6) 

Then, equations (I2.5P can be rewritten in the form 
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where B lj2 = (B ± 3cnF). Thus, the original problem (12. 3p is reduced to the problem 
of two independent pseudo-particles with the 'angular momenta' J\ and J 2 placed in 
the separate effective magnetic fields B\ and B 2 , i.e. the vectors J\ and J 2 uniformly 
rotate around the axes Bi and B 2 with frequencies u)\ = |Bi|/2c and u 2 = \B 2 \/2c, 
respectively. The quantization of subsystems like that is quite simple. In this case, 
the integrals of motion are the 'angular momenta' and projections of the 'angular 
momenta' Ji j2 onto the axes B lj2 , respectively. In the 3D case, the quantization of 
angular momenta gives (see also quantization condition f ll.4jl ) 

Jl = (j + 1/2) 2 , (2.8) 

where j is the 'angular' quantum number. Comparing the right-hand sides of equations 
(12.61) and (12. 8ft the following value of angular quantum number is obtained j = (n — 1)/2. 
For quantization of projections the well-known result is valicfl 

^(Ji • Bi) = ntUt, — (J 2 ■ B 2 ) = n 2 u 2 , n 12 = -j, -j + 1, j - 1, j. (2.9) 
2c 2c 

The first correction to the energy E\ is the perturbation XV (r, v) = — (F ■ r) + (B ■ L) 
averaged over the Kepler period. Employing (I2.9p . it can be written in the form 

Ei = Uiiii + uj 2 n 2 . (2.10) 

This result coincides with the quantum first-order correction [7], [8]. 

In the case of crossed electric and magnetic fields uj% = u 2 . The first-order correction 
is degenerated like in the case of Stark or Zeeman effects and an individual state cannot 
be defined. This degeneracy is removed in the second order of perturbation theory. This 
theory was developed in quantum approach only j9]. However, a classical approach can 
be developed in the same manner as in the case of the quadratic Zeeman effect, which 
is presented in the next section. 

2.1.2. Quadratic Zeeman effect The problem of a hydrogen atom in a magnetic field 
has fundamental importance. Attention to this problem significantly increased after the 
discovery that energy splitting at avoided crossing between adjacent manifolds decreases 
exponentially when the principal quantum number n increases |10j . 

In contrast to a hydrogen atom in an electric field the problem with a magnetic 
field is not separable. For the sake of definiteness, we choose the orientation of the 
homogeneous magnetic field B along the z-axis. Since the Hamiltonian of the hydrogen 
atom in a magnetic field (p 2 = x 2 + y 2 , u = B /2c is the cyclotron frequency) 

H = p^L + ^f + L , L , (2 . n) 

Because of non-complete understanding of the quantization rules that time, in [2], [5] and [6] the 
erroneous value of the angular quantum number j was ascribed as j = n/2, i.e. the quantum numbers 
ni t 2 took the semi-integer values instead of integer and vice versa. 
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Figure 1. Effective potential i7 e //($) as a function of i? for A = —0.5, /i = 0.1 (a), 
and A = 2, ^ = (b). 



is invariant under rotation around the 2-axis, the L z -component is conserved. Then, the 
motion along the azimuthal angle tp is separated out and the semiclassical quantization 
condition along (p gives 

L z =m, (2.12) 

where m = 0, ±1, ±2, ... is the magnetic quantum number. After the separation of the 
azimuthal angle ip the problem is reduced to a 2D non-separable one. 

At B = the electron moves on the Kepler elliptic orbit. In this case, the angular 
momentum L and the Runge-Lenz vector A are the two additional (to the energy) 
integrals of motion. Under action of the weak magnetic field these integrals start to 
change in time slowly, so that in the first order of perturbation theory the combination 
($ is the angle between the Runge-Lenz vector A and the z-axis) 

A = 4A 2 -5A 2 Z = A 2 (4-5cos 2 tf), (2.13) 

is conserved Taking into account that < A 2 < 1, the range of A values is 

— 1 < A 2 < 4. For A = the Runge-Lenz vector lies on the surface of the double cone 
fl specified by the condition cot^o = 2 (see Fig.l). For A ^ all trajectories split 
up into two classes: the trajectories with A < librate inside the double cone Q and 
trajectories with A > librate outside this cone. Thus, all states are localized in two 
non-overlapping domains. This unique property leads to the effect of the exponential 
smallness of energy splitting at avoided crossings between adjacent manifolds discovered 
by Zimmerman et al [TO] . 

The generalized momentum conjugated to the coordinate d is the angular 
momentum component perpendicular to the plane passing through the z-axis and the 
Runge-Lenz vector A (fi = m/n) 



=nJl+ A - J* = ny/l-U eff (0). (2.14) 
V 1 — 5 sin 1/ sin v v 
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Fig.l shows the effective potential U e ff(d) for two different cases: A < and A > 0. 

For negative values of A, the motion in intervals [$i,$2] and [$3, $4] is classically 
allowed (see Fig. la), the value of A is double degenerate and the final expression for the 
quantization rules has the form [T2] 



1?2 

Lj_(#)d0 = 7r(s + l/2), 

L L (-d)M = vr(s + 1/2), s = 0,1,2,... (2.15) 

Equations (12. 15ft give two degenerate states which are symmetric and antisymmetric 
to the xy plane. In the case of m = the turning points *&i and $4 do not occur. 
Instead, the two singularities appear at 1? = 0, tt due to the caustic at the z-axis which 
is characterized by the same Morse index as an ordinary turning point (see, e.g., [38] 
§49). Therefore, the quantization conditions in this case are formally obtained from the 
quantization conditions (12. 15ft by setting d\ = and $4 = n. The analysis of the roots 
of the function L±("&) shows that the states with A < exist only when m < n/yh. 

When A > 0, the region of a classically allowed motion is the interval [$5, $6] ( see 
Fig. lb). The A values in this case are nondegenerate and are determined from the 
quantization condition [12] 



[ * L±0&)d'& = n(k + 1/2), fc = 0,l,2,... (2.16) 

These states are localized outside the double cone Q and their parity with respect to 
the xy plane is equal to (— l) fc . 

In the first (with respect to u> 2 ) order of perturbation theory the quadratic Zeeman 
energy shifts are expressed in terms of the scaled value p 2 /n 2 averaged over one period 
of the Kepler orbit 



e = ^^= 1 -(l + ^ + A). (2.17) 

The A values are determined by the quantization conditions (I2.15P and (I2.16p . Of 
greatest interest are the outmost levels in a given {nm}-manifold, since they are the 
first to undergo overlap in the course of the approach to each other of two neighboring 
manifolds with increasing strength of the magnetic field. These energy levels correspond 
to the lowest levels in the effective potential U e ff($) in the quantization conditions (I2.15P 
and (I2.16p . For these states the potential U e ff($) can be approximated by the harmonic 
oscillator and the scaled energy shift inside the cone Q is (A < 0) 

e osc . = aV5 + 25a 2 + y/Efi - 5a 2 (2. 18) 

and outside the cone ft (A > 0) 



Base. = \~^+^-^ + K 2 - ^fi 2 . (2.19) 

where a = (2s + ^)/n and k = (2k + l)/n. In Table 1 the scaled energy shifts (I2.18P and 
(I2.19P are compared with the quantum calculations. The quantum results have been 
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Table 1. Comparison of the quantum results e q with the semiclassical approximation 



e osc . (Eqs.dHTU) and (EUg) ). 



n = 40, m = 



n = 40, m = 4 





0.055 0.053 2.45 2.44 

1 0.159 0.142 4 2.03 2.06 

2 0.255 0.212 8 1.65 1.75 

3 0.342 0.267 12 1.32 1.51 



0.251 0.276 2.43 2.43 

1 0.338 0.365 2 2.22 2.22 

2 0.415 0.436 4 2.01 2.04 

3 - - 6 1.82 1.88 



obtained by diagonalisation of the energy matrix within the given {nmj-subspace. The 
agreement for the position of levels is very good for low-lying states. For higher states 
the agreement becomes less satisfactory since the applicability of harmonic oscillator 
approximation breaks down. 

With the increase of the magnetic field the first avoided crossing arises between 
the lowest-energy state |1) of the n + 1 manifold and the highest-energy state |2) of the 
n manifold. In the classical approach instead of avoided crossing we obtain the exact 
crossing of energy curves since the first (A ~ — 1) and the second (A ~ 4) state are 
located in two nonoverlapping regions of configuration space (see Fig.l). The splitting 
is obtained in the quantum approach as the matrix element of diamagnetic potential 
between the first |1) and the second |2) state: 



In the first order of quantum perturbation theory the Schrodinger equation is separable 
in elliptical-cylindrical coordinates on a sphere in a four-dimensional momentum space 
[TTj [T2"| IT3] . Using uniform semiclassical approximation for the wave functions of the 
states |1) and |2) in this coordinate system (see, e.g., [E]), the splitting is obtained in 
the form [12] 



which coincides with the probability of the under-barrier penetration in the effective 
potential f/ e //(^) with A = — 1 in the interval $ G [^2,^0] ( see Fig-la); an d with 
A = 4 in the interval $ e [^0,^5] (see Fig. lb). Figure 2 demonstrates perfect agreement 
between this result and experimental data [TO] . 

2.1.3. Helium atom; equivalent electrons The first attempts to develop a semiclassical 
perturbation theory for a helium atom were made in the old Bohr theory (see, e.g., 
review [IE]). According to heuristic concepts accepted at that time, only the simplest 
symmetric trajectories were considered, which is wrong from the modern point of view . 
The self-consistent perturbation theory was developed for equivalent electrons (having 
the same principal quantum number n) with total angular momentum equal to zero [17] . 



AE = (1|<V|2). 



(2.20) 



AE ~ exp{-nln[(v / 5 + 2)(v / 5+ l)/2]} ~ exp(-1.92n) 



(2.21) 
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Figure 2. Energy splitting between the lowest-energy state of the n + 1 manifold and 
the highest-energy state of the n manifold as a function of principal quantum number 
n: solid line - the approximation (|2.21j) . solid squared - experimental data from [10] , 



A distinguishing feature of the classical perturbation theory for equivalent electrons 
of a heliumlike system with the nuclear charge Z is the presence of the accidental 
degeneracy in the unperturbed stat^B In this case, the perturbation series is the 
series over half-integer power of the small parameter A = 1/Z (see [18] §108). To 
construct this series, first of all, the proper variables should be introduced. When the 
total angular momentum is equal to zero, the trajectories of the two electrons are on 
the same plane, regardless of the magnitude of the electron-electron interaction, and 
their angular momenta are equal but oppositely directed. In the zeroth order (without 
electron-electron interaction) both electrons move along the Kepler ellipses. The mutual 
orientation of the ellipses in the plane is specified by an angle $ = $i — d 2 with t?j being 
the angle between the Runge-Lenze vector of the i-th electron and z-axis (see Fig. 3). 
The position of the i-th electron on the ellipse is determined by the Kepler anomaly 
which is related to time by the relation 

t-n = - ^l + 2E i LyZHmi l ) (2.22) 

where Xj is the instant of passage through the perihelion, Ei is the energy, Li is the 
angular momentum and Tj is the period of the z-th electron. Since the dependence 
on time of the variables £i and £2 is not independent, one of them (for the sake of 
definiteness £ = £1) and the time delay r = t\ — r 2 between the passage of the first and 
second electrons through the perihelion should be used. 

+ Accidental degeneracy means the commensurability of the oscillation periods for two or several 
coordinates, which takes place not always but at some initial conditions. 
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Figure 3. Trajectories of two equivalent electrons for the case Xm = n/lO. From |17j . 

In the planar case (L = L\ — —L2), the motion of both electrons is described by 
eight Hamiltonian equations. It follows from four of them that the quantities (Li — L 2 ) 
and (■#! + $2) are exactly conserved, while (Ei + E 2 ) and (ri + r 2 ) are conserved in the 
first-order. A nontrivial role is played by the equations 

where Eq = —Z 2 /2n is the unperturbed energy of one electron, AE = (E% — E 2 )/2 and 
V(l;, L, Ei,E 2 , t) is the electron-electron Coulomb interaction. The small parameter A 
can be eliminated from equations f l2.23|) and f l2.24[) by changing from AE to e — AE j \f\ 
and introducing new 'times' s\ = y/\t and s 2 = Xt for the first and second pairs of the 
equations, respectively. It follows hence that the rate of change of r and e is of the 
order of V% and that of L and d of the order of A, i.e., the angular momentum L and 
mutual angle •& vary infinitely slowly compared to time delay r and scaled energy e as 
A — > 0. In addition, the energy transfer AE is always small of the order of v^A, therefore 
the change of the parameters E\ and E 2 , as well as the periods in the electron-electron 
interaction should therefore be neglected, and these arguments of V will hereafter be 
omitted. Since the characteristic frequencies are different, the motion along variables 
e, t and L, -d is adiabatically separated. 

Most sensitive to electron-electron interaction are the variables e and r. They 
oscillate with period of the order of A" 1 / 2 , during which L and i? can be regarded as 
constant. On the other hand, the oscillation over the Kepler anomaly £ has a high 
frequency compared to e and r. Replacing the interaction V by its averaged value over 
common period (T = T\ = T 2 ) 

V(r; L,&) = ^ j T V(C{tl L, 0, r)dt (2.25) 
1 Jo 

we obtain one-dimensional problem in which L and enter as parameters. The 
important feature of this one-dimensional problem is that only the ground state of 
the potential AV(r; L, 0) corresponds to the unperturbed state and for quantization the 
corresponding Schrodinger equation should be used. Since AE and r are canonically 
conjugated variables, the Schrodinger equation reads 

3 "' +\V(t;L,#)U(t',L,#) =£(L,t)iP(t;L,0) (2.26) 



2E dr 2 
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Figure 4. Effective Hamiltonian as a function of v and \. The dashed line is the 
function L(x) for Xm — 27r/5 and q = 0.61. From [17] . 



with the periodic boundary condition: ip(0; L, i?) = ip(T; L, •&). For all excited states of 
equation (I2.26P AE is finite at A = and they do not go over into the states of the 
unperturbed problem. In the first order of perturbation theory the wave function of the 
ground state is constant, which is ipo — 1/VT. The first-order correction 

£(L, 0) = f T ^(r)V(r; L, $)dr = - f V(r; L, $)dr (2.27) 
J o i jo 

is additional integral of motion which is an effective Hamiltonian in the angle After 
the substitution of (12. 25ft into (12. 2 7ft this expression takes the form 

S(L,,)='f T f ] * l *> t ,. , (2.28, 
T Jo Jo |ri(ti) - r 2 {t 2 )\ 

which coincides with classical averaging but for the nondegenerate case. This leads to 

the important conclusion that from the point of view of the quantization conditions 

there is no difference between the nondegenerate and accidental degenerate motion. 

Using the explicit expression for the Kepler ellipses the scaling with respect to the 

principal quantum number n is obtained 

8(L,#) = ^v(u,&) (2.29) 

where v = L/n. The numerically calculated effective Hamiltonian v(x,v) is shown in 
Fig.4. 

The dependence of L on -d is determined by the condition that £ (L, #) is a constant 
on the trajectory, i.e. 

v(u,d)=w (2.30) 

The constant w is determined from the quantization condition, which (in terms of the 
variable \ = (n — $)/2) is the same as considered for a hydrogen atom in a magnetic 
field at m = (see Sec. 2. 1.1). Finally, the quantization condition takes the scaled form 

Xm ijr 

Hx)dx = ^ ( 2 - 31 ) 
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ff,V 0,8 % 

Figure 5. Scaled correction to the unperturbed energy: solid line - classical 
perturbation theory, open circles - quantum perturbation theory, solid triangles - exact 
values for He, solid circles - exact values for B 3+ . The parentheses contain the state 
quantum numbers (n, k). From |17j . 

where Xm is the turning point, q = (2k + l)/2n and A; is a new quantum number 
(k — 0, 1, 2, n — 1). In the first order, the correction E^ 1 ' to the unperturbed energy 
is equal to the average electron-electron interaction 

= ^w(q). (2.32) 

Figure 5 shows a plot of w(q) and also the corrections, recalculated in accordance with 
the scaling rule f )2.32p . to the unperturbed energy in first-order quantum perturbation 
theory for helium and obtained by exact calculation for He and B 3+ [17]. As it is seen, 
the agreement improves with increasing nuclear charge Z. The discrepancy between 
classical and quantum perturbation theory becomes noticeable only at q > 0.5. It is 
attributable to the fact that as q — > 1 the turning point \ m approaches the singular 
point x — 7r /2 (see Fig.4). In this case, the quantization condition (I2.3ip must be 
replaced by the modified quantization condition in which simultaneous account is taken 
of the turning point and the singular point at \ = 7r /2- 

2.2. Adiabatic switching method; hydrogen atom in external fields 

The extension of the semiclassical quantization to the non-separable systems was 
proposed by Einstein in 1917 [2T] . who pointed out that the classical trajectory is 
Lagrangian manifold, so if the classical trajectory forms a torus in phase space, then 
in the N dimensional problem only N topologically independent quantization contours 
on this torus exist, and the actions along these contours are, according to the Liouville 
theorem, invariant with respect to the deformation of the contours. However, for a 
system with several degrees of freedom the calculation requires an excessive amount 
of computer time, mainly spent in rejecting unsuitable trajectories. To avoid this 
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problem, the adiabatic switching method was proposed [20] • This method is based 
on the quantum Born-Fock adiabatic theorem [22] and the correspondence principle 
between quantum and classical mechanics^]. The adiabatic switching procedure is quite 
simple. To calculate energy spectrum of non-separable Hamiltonian H, first of all, we 
need to guess a reference Hamiltonian Hq which is, on the one hand, solvable (separable) 
and, on the other hand, has the same topology of caustic as Hamiltonian H. Then one 
have to compute with the classical equations of motion the development in time of the 
initially quantized trajectory during a slow switching on the interaction V = H — Hq 
(if (A) = H + \(t)V). When the interaction has been fully switched on (A = 1), the 
quantized trajectory for Hamiltonian H and corresponding eigenvalue of energy are 
obtained, as well as, for all intermediate values of switching parameter A. It is expected 
that the more slowly the interaction is switched on, the more precisely the quantization 
conditions are satisfied. Of course it is heuristic technics without proof but it works in 
many cases; it is a simple and effective tool to calculate the energy spectrum of involved 
mult i- dimensional systems. The only thing necessary to control is that the topology of 
caustics does not change during the switching [20] . 

In the case of hydrogen atom in crossed electric F and magnetic B fields, the 
reference Hamiltonian Hq is the Hamiltonian of the hydrogen atom [25]. Due to the 
high dynamical symmetry of the Coulomb interaction all trajectories in configuration 
space are closed lines - ellipses and just one quantization condition can be written which 
specifies principal quantum number n only. When external perturbation is introduced, 
the trajectory fills up 3D domain in the configuration space and missed quantization 
conditions can be formulated. The perturbation theory for hydrogen atom in crossed 
electric and magnetic field is presented in Sec. 2. 1.1. Thus, the initial conditions for 
the quantized trajectory is determined by conditions (12.91) . Next step is numerical 
calculation of classical equations for electron - 

d 2 r r wxi",, 1/ ld\(t) . 

— + - = -A t [F + -(v x B + -—^-A (2.33) 
dt z r 6 1 c J c at 

- during the switching external fields from A = to the final values A = 1. In equation 

( 12.331) A is the vector-potential (B = rot A). The last term on right-hand side of 

equation (I2.33P is the additional force arising as a consequence of the time dependence 

of the magnetic field. Although this force, being proportional to the switching rate, 

disappears in the adiabatic limit, it cannot be neglected. It ensures the adiabatic 

invariance, e.g., in the case of a pure magnetic field (F = 0). Obviously, if it were 

neglected, the electron energy would be conserved exactly, that is wrong. 

Figure 6a shows the energy levels (full curves) of the ground state as a function 

of field strengths, for different combinations of parameters a = B/cF and the angle 

9 between F and B. As a limit, quantum- mechanical results (chain curves) for the 

ground-state energies of a hydrogen atom in a magnetic field only (E = 0) and in 

f Often, in the literature this method is erroneously associated with the Ehrenfest adiabatic principle 
|23) . The Ehrenfest principle was formulated for separable problems only. Moreover, it was proved that 
in the non-separable case it breaks down (see, e.g., [21]). 
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a) b) 

Figure 6. a) Ground state energy levels as functions of field strength for various 
combinations of parameters a — B/cE and cos 0. Chain curves are the quantum results 
for the cases of pure electric [27J and magnetic [53] fields. Broken curves demonstrate 
the non-adiabatic evolutions of the system, b) Energy levels of excited states, 
(n, ni,n 2 ) = (3,-1,-1) (the levels shifted downwards) and (n,ni,n 2 ) = (3,1,1) 
(the levels shifted upwards), as functions of field strength with constant parameters 
a = 9, 27 and cos# = 0.5. Broken lines are the first-order perturbation theory results 
(|2.10p . The arrows indicate ionization in both figures. From |25) . 

an electric field only (B = 0) are shown. The curve labeled by E = is the result 
|26j, obtained by diagonalising the energy matrix in an extensive basis of hydrogen 
wavefunctions. The curve labeled by B = is a quantum fourth-order perturbation 
theory result (see, for example, [2?]). Figure 6b shows the results (full curves) for the 
states defined by the weak-field quantum numbers (n,ni,ri2) = (3,-1,-1) (the levels 
shifted downwards) and (n,rii,n 2 ) = (3,1,1) (the levels shifted upwards). The angle 
between the fields has been kept constant (cos^ = 0.5) and the results are presented 
for two different ratios of field strengths: a = 9 and a = 27. Unlike the ground-state 
case, here the linear shift, given by ( 12 . 1 0[) (broken lines), dominates in the limit of weak 
fields. The state (3,-1,-1) is bound stronger and the ionization limit was not reached in 
the range of fields shown in figure 6b. For the (3,l,l)-state the non-adiabatic behaviour 
was found only very close to the ionization limits. In Fig. 6 electric field scale is common 
to all curves, while the magnetic field scale corresponds to curves labeled by a = 9 for 
the ground state and a = 27 for excited states. In all cases, the switching rate used in 
calculations was d\(t)/dt = 2.5 x 10~ 5 au. 

Again, adiabatic switching method has certain advantages with respect to the 
straightforward semiclassical quantization of non-integrable systems. It is free of such 
non-local problems as finding the caustics of the classical system and searching for the 
initial conditions which correspond to quantized trajectories (see, for example, [28J ) . In 
addition, this method provides all intermediate energy values at < A < 1. 
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Table 2. Classical lifetime t c \ (Eq. (|2.35[) ) and quantum lifetime T q [34] for spontaneous 
decay of excited states of a hydrogen atom (in nanoseconds). 



nl 


2p 


3p 


Ap 


5p 


6p 


3d 


Ad 


5d 


6d 


4/ 


5/ 


6/ 


T C 1 


1.68 


5.66 


13.4 


26.2 


45.3 


15.7 


37.3 


72.8 


126 


73.1 


143 


247 


T 1 


1.60 


5.27 


12.3 


22.2 


40.8 


15.5 


36.2 


69.7 


119 


72.5 


140 


240 



3. Spontaneous decay of excited states of a hydrogen atom 

A very interesting situation takes place in the problem of spontaneous decay of excited 
states of a hydrogen atom. From the classical point of view this process is due to the 
bremsstrahlung (see, e.g., [32j [33] ) . In the old Bohr quantum theory the bremsstrahlung 
was the main obstacle to the description of the stable atomic ground state within the 
framework of classical mechanics since the accelerating electron should emit radiation 
losing energy and finally should fall down onto the nucleus. However, it does not 
contradict the evolution of the excited states which are unstable. 

According to classical electrodynamics, the rates of decrease of the energy E(t) and 
the angular momentum L(t) are described by the system of equations (see, e.g., [35] 
§75) 



dE(t) {-2E(t)}y\ o , . 



dt ' 3c 3 L{t) 
dL(t) 2{-2E(t)Y' 2 



{3 + 2E(t)L 2 (t)}, 

(2.34) 



dt 3c 3 L 2 (t) ' 

The lifetime r is defined by the expression P = const x exp(— t/r), where P is the 
population of the excited state. In classical mechanics there is no such concept, but we 
can estimate it as the time t c \ in which the initial angular momentum decreases by unity 
(AL = 1), which corresponds to the emission of one photon. In the leading order of 
classical approach (n, I — > oo) we can neglect the dependence on time in the right-hand 
side of equations ( 12.34)) . Then, after the substitution for the energy E = —l/2n 2 and 
angular momentum L = (I + 1/2) we obtain the following estimate of lifetime in terms 
of the quantum numbers n, I [33] 

Td = ^c 3 n 3 (l + 1/2) 2 = 9.32 x 10" n n 3 (/ + 1/2) 2 sec (2.35) 

In Table 2 these results for several states are presented in comparison with quantum 
calculations of total lifetimes r q of decay from the nZ-state into all low-lying states [34J . 
This table demonstrates amazing agreement even for the low values of quantum numbers 
n and /, though the lifetime is not well-defined in classical theory. 

The system of equations (12. 34ft has an exact analytical solution but, according to the 
correspondence principle, only the leading classical order with respect to n, I — > oo has 
physical meaning. The next terms compete with unknown quantum corrections. Since 



CONTENTS 



17 



the corrections of the same order caused by the different sources usually compensate 
each other, we obtain worse results taking them partly into account. It is the reason of 
bad agreement obtained in [32] where classical and quantum orders were mixed up. 

Classical physics is deterministic, that is why the decay process has one channel 
providing the total lifetime only, and there is no way to analyze state-selective lifetime. 

Generally, decay processes are not well-defined in quantum mechanics. Here we 
can not separate the quantum system from measurement. From the very beginning, 
it is formulated as a exponential decreasing in time of the population of state which 
is rather a classical treating of the problem. It is not clear how to determine the 
width of the energy level in the degenerated case, e.g., in the case of hydrogen atom in 
excited state at n — 2. If we add weak spherical perturbation, the degeneracy is moved 
and correct wave functions are the spherical functions \2s) and \2p). Since the dipole 
matrix element between \2s) and ground state | Is) is equal to zero, \2s) is a meta-stable 
state having macroscopic lifetime - 1.4 sec. The second state \2p) has typical atomic 
lifetime 1.6 x 10~ 9 sec. If we add weak electric field the correct wave functions are 
characterized by parabolic quantum numbers n 1; n 2 , m (n = n\ + n 2 + m + 1). At 
n = 2 and m = there are two states \nri1n2) = |210) and 1 201 ) which have the same 
lifetime of order 2.3 x 10 -9 sec. But what is the lifetime for the hydrogen without any 
perturbation? There are huge difference between lifetime in spherical and parabolic 
coordinates. Here we have an unsolvable conflict with the superposition principle which 
is typical of classical processes. 



2.4- Gutzwiller's approach; broad resonances 

The contribution of the unstable periodic orbit to the trace of the Green function is 
determined by the Gutzwiller formula [52] 

iT{E) - exp{zn[S(E)/h-\ir/2]} 

9{E) ^~^Tti sinhK(£)/2] (2 ' 36) 

where S(E), w(E), T(E) and A are the action, the instability exponent, the period 
and the number of focal points during one period, respectively. After the expansion 
of the denominator, according to [sinh(x)]" 1 = 2e~ x Y^kLo e~ 2kx , and summation of the 
geometric series over n one can see that the response function f!2.36j) has poles at the 
complex energies E ks whenever 

S{E ks ) = h\-n/2 - ihw(E ks ) ( k + ~) + 2snh, k, s = 0,1,2, ... (2.37) 



2, 

Equation (I2.37P is a transcendental equation with respect to E ks . This approach is not 
well-defined [T], but, in the case of the shortest period associated with unstable periodic 
orbit in separable system, the poles of (12. 36ft can be interpreted as manifestation of the 
resonances in the continuum. The width of these resonances is the linear function of h 
in spite of the under-barrier resonances whose width is exponentially small with respect 
to h — > ( oc e~ 2 ! s !/ n , where S is under-barrier action). The Gutzwiller's approach was 
applied to the realistic systems - the scattering of electron on Coulomb potential in the 
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Table 3. Real parts Eq and widths Tq (in atomic units) for a sequence of resonances 
with parabolic quantum numbers (n\ , n% , m) at field strength F = 8 kV/cm, as 
obtained from Gutzwiller's approach (|2.44j) . E q and T q are the results of the full 
quantum calculations |58j . 

(nx,n 2 ,m) (23,0,0) (23,1,0) (23,0,1) (24,0,0) (24,1,0) (24,0,1) 

£ G x 10 4 1.899 1.774 2.597 3.360 3.251 4.021 

E q x 10 4 1.949 2.039 2.698 3.382 3.433 4.090 

r G x 10 4 0.499 1.480 1.145 0.632 1.887 1.355 

T q x 10 4 0.524 1.681 1.188 0.638 1.955 1.369 



presence of magnetic [531 IS] an d electric fields [551 [56], an d scattering on two- Coulomb- 
center potential [56] . 

As an example let us consider resonancees in elastic scattering of electron on the 
Coulomb centre in the uniform electric field F [56]. The Hamiltonian is separable in 
parabolic coordinates u, v and if (x 2 + y 2 = u 2 v 2 , 2z = u 2 — v 2 and u, v > 0). After 
introducing a new time-like variable r defined by the relation dr = dt/(u 2 + v 2 ), the 
problem is reduced to two separated one-dimensional systems: 



2u 2 



H i = ^+ '^6r ~ En 1 + ^Fu 4 = 2/3, (2.38) 
H 2 = ^ + - Ev 2 - ^Fv A = 2(1 - 0), (2.39) 



where /3 is a separation constant. At positive energy E the motion in the coordinate u 
is bounded whereas in the coordinate v it is unbounded. The system f l2.38p . fl2.39p has 
hyperbolic fixed point whose location in the v — p v space at energy E = Ek s is given 
(up to first order in Ti) by 

J^\h(2E ks )- 1 / A e l7T/ \ p* v = 0, 0* = l + i\m\hJE ks /2 (2.40) 



This fixed point corresponds to unstable state since in its vicinity the Hamiltonian H2 
has a form of turned oscillator (8p v = p v — p*, Sv — v — v*) 

5 2 H 2 = ^(5p v ) 2 -AE ks (Sv) 2 - (2.41) 

In the first order with respect to % — > the corresponding action and instability exponent 
are 

r u o 1 i 

S(E ks ) = 2 / j2E ks u 2 + 4 - FuHu - \m\ith - -\m\hw(E ks ), (2.42) 
Jo 2 

w{E ks ) = ApE^J =, (2.43) 

v Jo y/2E ks u 2 + 4 - Fu 4 

where «o is the external turning point. The substitution ( 12. 42ft . (12.431) into ( 12. 3 7ft gives 
the transcendental equation for complex energy E niTl2 (s = ni, k — n 2 , A = 2) 

S{E nin2 ) = (2n x + l)Trh + ih(n 2 + -)w(E nm2 ). (2.44) 
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which is, up to the first order of h, equivalent with equations (15)-(17) in Ref.joT] 
obtained in quantum approach by the comparison equation method. Table 3 
demonstrates a good agreement between Gutzwiller's approach ( 12.44ft and full quantum 
calculations [58] for the series of resonances which lie close to the real E-axis. They 
all correspond to large values of parabolic quantum number n\ and small values of ri2 
and m. The widths Tq are comparable with the differences of the positions Eq of the 
adjacent resonances. In this sense they can be referred to as "broad resonances". They 
explain, e.g., the experimentally observed resonant structure of photoionization cross 
section of the hydrogen atom near the ionization limit in strong electric fields [581 EH] • 

This problem was considered by Wintgen too [55]. His result, although qualitatively 
good, does not give correct values for the complex poles of expression f!2.36|) in the case 
m = 0. The approach presented in [56] is based on the complex periodic orbits in the 
m^O case and gives the correct result for m = as m — > limit. 

The same approach with the same result was developed for broad resonances in the 
scattering of electron on two-Coulomb-center potential [56] . 

2. 5. Poincare section; resonances in helium atom 

The Poincare section is a powerful tool to study the irregular motion in 2D non-separable 
problems because of the visualization of a dynamical system on one plot. In this case, the 
phase space is 4D with two coordinates (x, y) and two momenta (X, Y). At fixed energy 
E one of these variables can be expressed from the equation H(x, y, X,Y) = E as a 
function of the energy E and the remained variables. Let it be, for instance, the second 
momentum Y(x,y,X,E). Now if we plot on the plane {x,y} the values of variables 
x(t) = Xi and X(t) = Xi each time ti when the coordinate y(t) takes some value yo 
then, in the general case, the distribution of points (xj,Xj) in this plane is irregular. 
But if an additional integral of motion (may be approximate one) exists - A(x, y, X, Y), 
then from the transcendental equation A(x,y , X,Y(x,y , X, E)) = A we obtain the 
functional dependence of x on X and constants yo, E and A: x(X; yo, E, A), where A is 
the value of the second integral of motion. As t — > oo (i — > oo) the points (xi,Xi) from 
one trajectory fulfill this line. Calculation of many trajectories gives the Poincare section 
of the system, which visualizes the information about regular and irregular regions of 
motion. In the regular region there can exist two types of fixed points, elliptic and 
hyperbolic, associated with stable and unstable periodic orbits, respectively. 

Consider a collinear arrangement of a nucleus Z and of two electrons, both being on 
the same side of the nucleus. It is a two-dimensional problem (position of inner and outer 
electrons on the axis), where we can use the Poincare section which is shown for helium 
in Figure 7a. In this section, the phase space position (r 1; pi) of the outer electron is 
monitored each time the inner electron approaches the nucleus (v2 = 0). From Figure 7 
one can see the elliptic fixed point at p\ = and r\ ~ 7, which corresponds to the two- 
electron periodic orbit (Fig. 7b upper trajectory). Such a motion is a coherent oscillation 
of both electrons but with large difference in their positions and velocities. The outer 
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Figure 7. a) The Poincare section for collinear configurations of the electrons. The 
coordinate r± and the momentum p\ of the outer electron is monitored whenever 
the inner electron approaches the nucleus (r% = 0). b) The periodic trajectory 
corresponding to the elliptic fixed point on the left panel (a) and a nonperiodic but 
regular trajectory in its neighbourhood (b). From [31] . 



electron appears to stay nearly frozen at large distance. This localization of outer 
electron is amazing since the total charge of the core (the nucleus and inner electron) is 
attractive. The existence of such states is explained by the resonant interchange of the 
energy between two electrons, i.e. it is a purely dynamical effect - at the conventional 
analysis we should expect that outer electron should fall down onto the nucleus since 
the total charge of nucleus and inner electron is equal to -1. Thus, it is very surprising 
that these classical configurations are extremely stable. 

Semiclassical quantization near the periodic orbit is based on the Gutzwiller's 
formula (|2.36|) . but, since the orbit is stable, the instability exponent must be replaced 
by the stability parameter 7 = — iw. The final expression of the quantized classical 
action of the tori surrounding the periodic orbit is given by [31] 

SnM = s + 1/2 + 2(^+1/2)7! + (I + l/2) 72) s,M = 0,l,2,... (2-45) 

where 71 = 0.46164 and 72 = 0.06765 the winding numbers, which describe the 
behaviour of nearby trajectories. The quantum numbers s, k and I reflect the 
approximate separability in the local coordinates parallel and perpendicular to the 
periodic orbit; s is the quantum number along the periodic orbit, k represents bending 
degree of freedom and I corresponds to the perpendicular degree but preserving 
collinearity. Using the scaler property S = S sc /J{ — 2E), the quantized energy levels 
are obtained 

S 2 

Enkl = ~ [s + 1/2 + 2(k + 1/2)71 + (/ + l/2) 72 ] 2 ' (2 ' 46) 
where S sc = 1.4915 is the (scaled) action of the periodic orbit for helium. 

Table 4 shows the semiclassical energies (I2.46P in comparison with full quantum 
calculation of energies E q . In the table the width of the states T q with respect to the 
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Table 4. Scmiclassical energies E sc i from expression ()2.46|) . quantum energies E q and 
widths T q /2 for the states of the series (s, 0, 1)). From [31] . 
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autoionization is also presented. One can see extremely small values of the widths for 
this series of asymmetric double excited states. 

3. Ion-atom and electron-atom collisions 

3.1. Electron-impact detachment of negative ion near the threshold 

The ionization of a negative ion by an electron impact draws much attention because of 
a huge cross section. This reaction is very important, for instance, in understanding of 
stellar atmospheres. 

First of all we need a classical description of the negative ion. The classical model for 
interaction of electron with neutral core of negative ion is based on the same physical 
assumptions as the zero-range potential model in the quantum theory [37], i.e. the 
potential well supporting one weakly bound s-state is spherical, narrow, and deep. The 
electron is oscillating along the diameter between the opposite turning points d and —d. 
The frequency of this oscillation u is the large parameter of the classical model. The 
bound s-state is represented by the ensemble of such trajectories uniformly distributed 
in all directions. 

Now let us consider the detachment of the negative ion by an electron as a projectile. 
Near the threshold, because of the Coulomb repulsion between a negative ion and a 
projectile, there is another large parameter - the distance of the closest approach of the 
projectile to the negative ion Rq. Since Rq is large, there is no interaction between the 
projectile and neutral core. Then, the equation of motion for the projectile is 

dt 2 I?"! — r 2 \ 6 

where r± is the position vector of the projectile and r 2 is the position vector of the 
bound electron. Since r 2 < d, r\ > Ro and Ro 3> d, we can expand the right-hand side 
of equation (13. ip over the multipole series. For our purposes, it is sufficient to consider 
the dipole approximation 

d 2 ri = n _ ra + 3(ri ■ r 2 )r 1 

dt 2 r\ r\ r\ 

The first term on the right-hand side of equation f 1 3 . 2 H describes the interaction with the 
fixed Coulomb center, but the second and third represent the rapidly oscillating force 
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with frequency u. Kapitsa developed an approximation method for this situation (see, 
e.g., [38], §30). In this method, the position vector of the projectile splits into two terms 
7*1 = R + r, where R is a smoothly-varying function and r is the rapidly oscillating 
part having zeroth mean value over a period 2u /ui. In the leading order with respect to 
u — > oo this method leads to the following results: 

l- 2 Eb 2 1 

2 R + W + R = E > 



1 



r 1 _ 3 (r 2 -R)R 



-B? R 3 ^ (3 ' 3) 

where b is the impact parameter and E is the initial energy of the projectile. Ionization 
occurs when the work done by the projectile A becomes equal to the binding energy of 
the negative ion e: 

A=[ ^^tp = e (3.4) 
J \r 1 — r 2 \ 6 

It is easily verified with the aid of f)3.3p that the work done by a projectile is an oscillating 
function whose mean value is zero, and the amplitude of the oscillation increases in 
the course of the collision. Consequently, the energy necessary for ionization is not 
accumulated gradually in the collision process, but is transferred in a quarter of the 
period while the bound electron moves from the nucleus in the direction opposite to 
that of the projectile. The orientation of each trajectory in the ensemble is characterized 
by spherical angles if with respect to the vector R(t). The classical trajectories, for 
which escape is possible, are then restricted to a cone Q(t) with •& < i? m (i). The size of 
the cone is obtained from (13.41) : 

rfcos^ m (t) _ 

R 2 (t) ~ [6 - b) 

The detailed analysis shows that the cone corresponding to the distance of the closest 
approach Q absorbs all previous cones Q(t). Therefore, the probability of ionization at 
a given impact parameter b is 

W(b) = —[ sin&d&dip = 1 - 4|r(l + Vl + AE 2 b 2 ) 2 (3.6) 
2tt Jfi AE 

where E t h = Je/d is the classical threshold of detachment. The normalization factor 
in equation (13 .6p is taken as l/2n because the opposite directions represent the same 
trajectory. The total detachment cross section is [39j 

a = 2nj o bm W(b)bdb = 2^i(£ - E th ) 2 (E 2 + 2 -EE th + ^E 2 h ), (3.7) 



where b m = JE(E — E t h)/(E t hE). The classical threshold energy E t h differs from the 
quantum threshold energy e, because, in the classical approach, the projectile needs extra 
energy to overcome the Coulomb repulsion and to reach the distance where ionization 
occurs. For e < E < E t h the detachment mechanism is purely quantum, it is an the 
underbarrier penetration. 
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Figure 8. The detachment cross section of H (or D~) as a function of electron 
impact energy. The full curve is the classical cross section (|3.7j) . the open circles are 
experimental data from [45] and the solid circles are experimental data from [46] . 

Figure 8 shows the cross section (13. 7\\ for electron- impact detachment of H (or D~) 
in comparison with experimental data. The parameters for H~ were taken from [37J: 
d = 2.7 and e = 0.0278. The classical results are in good agreement with experimental 
data [15] and follow the shape of experimental data [16] in the threshold region. However, 
they are shifted approximately by 2 eV with respect to the latter. This shift is equal to 
the difference between classical (2.76 eV) and quantum (0.76 eV) thresholds. 

It is not clear how to reproduce the classical mechanism of detachment in quantum 
mechanics. The detachment cross section in the threshold region depends critically 
on the interaction between the projectile and the rapidly oscillating bound electron. 
In quantum mechanics, it is difficult to take this interaction into account rigorously, 
whereas in classical approach this can be done simply and naturally. In [3Qj the authors 
tried to modify the classical approach. However, the proposed "corrections" are beyond 
the selfconsistent asymptotic approach. For instance, it was taken into account that 
not all electrons, satisfying the condition d < $ m , have time to leave the negative ion 
during collision, and, at the same time, the under barrier penetration is considered, in 
fact, as an instant process. 

3.2. Non-adiabatic transitions via hidden crossings 

In slow atomic collisions the inelastic transitions between electronic adiabatic states 
occur at the place of the closest approach of adiabatic potential curves - avoided 



CONTENTS 24 

crossings^. Initially in the application of the adiabatic approach the narrow avoided 
crossings due to the under-barrier resonant interaction of two electronic states located 
on different nuclei were used. Later [12] the new mechanism of nonadiabatic transitions 
via the so-called 'hidden crossings' were discovered. Hidden crossings happen when 
an electronic energy level touches the top of an effective potential. In the classical 
description the full dimensional electronic trajectory collapses into an unstable periodic 
orbit at this place. The transitions via hidden crossings are dominant in atomic 
collisions. They provide a complete description of inelastic processes in atomic collisions. 
The classical/semiclassical theory of hidden crossings was developed in [35]. In detail it 
is discussed in the review paper 



3.3. Binary encounter approximation 

For inelastic transitions of atomic electron in the charged projectile-neutral atom 
collisions, the binary encounter approximation can be applied. This approximation is 
based on the assumption that inelastic cross sections can be obtained from the two-body 
(projectile-atomic electron) elastic cross section after averaging over the momentum 
distribution of atomic electron in the initial atomic state. This scheme can be used in 
both quantum and classical approach. 

Initially, the classical binary encounter approximation was developed for fast 
collisions (see, e.g., review paper [UJ and references therein). In [19], the restriction on 
impact velocity v p » v e , where v p and v e are the initial projectile and atomic electron 
velocities, was removed. The most interesting case is the case when the projectile is an 
electron. In this case, the electro-electron velocity transfer from the projectile to the 
atomic electron is given by 

2(v 3 b-2v) , , 

AVe = ' (3 ' 8) 

v 4 o + 4 

where v = v p — v e is the electron-electron impact velocity and 6 is a vector of impact 
parameter (vA.b). Then the kinetic energy transfer is 

\jpfu\ 2v P v eV 2 b sintf cos x + 2f 2 + 4(« ■ v e ) 

AE{b) = VvTl ' (3 ' 9) 

where # is the angle between v p and v e , x is t ne angle between vector of b and the 
plane in which vectors v p and v e lie. The boundary of the region of b values which lead 
to an energy transfer greater than e is specified by the condition AE(b) = e, and the 
differential cross section with respect to velocities of projectile v p and atomic electron 
v e for the transfer of an energy AE(b) > e is 

E{v p , v e , e) = j e(AE(b) - e)£b = - v 2 p v 2 e cos 2 0), (3.10) 

| According to the von Neumann- Wigner theorem [41], the exact crossing of two adiabatic potential 
curves of the same symmetry is forbidden. 
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v 2 + 2e are 



where Q(x) is the step- wise Heaviside function, and v' p = Jri* — 2e, v' e = 
the projectile and atomic electron velocities in the final state. Assuming the isotropic 
distribution over v e in the initial atomic state the averaged cross section is obtained as 



Vp + Ve ^ (Vp + Ve)' 

v p v e 



(Vp - v e 



( / 2 
2 P 



Vp ~ v 2 ) 2 



- 2 



+ v 2 + 2v'v' 



-4 



V p > V' £ 



V n < V' 



(3.11) 



V'pV' e ~ (v 2 + V 2 - 2v' p v' e j 

The final expression for the cross section a(v p , e) is found by taking an average of (13.111) 
over the distribution of electron in the atom q(r). For the spherical symmetric atomic 
potential V(r) the distribution in the nZ-subshell has the form jH] 

. . 1 



Qm(r) 



2nr 2 T n ip n i(r) 



(3.12) 



where p n i(f) = y2[E n i — V(r) — L 2 /2r 2 ] is the radial momentum of atomic electron in 
the initial state E n i, L = h(l + 1/2) and T n i is the period of oscillation over r. Then the 
cross section is expressed as an integral between the turning points r\$ 

2 r 2 „ , rzTz: , dr 

T n l Jri 

For the fast collision (v p — > oo), the cross section does not depend on the mass of 
projectile and has a general asymptote 

2nZ 2 



2 [ r z I 

cr n i(v P , e) = — / S(v p , \J2[E nl - V(r)],e] 



Pni{r) 



(3.13) 



a n i(v pi e) 



3v 2 e 2 



[3e + 2{E nl - V n i)], 



(3.14) 



where V n i is the average potential energy (for Coulomb potential V n i = 2E nl ) and Z is 
the charge of the projectile. At e — \E n i\ expressions (13. 131) . (13. 141) give the ionization 
cross section cr ion (v p ). 

To calculate the cross section for excitation to a state with a principal quantum 
number n', e in (I3.12p should be presented in the form e = E n > — E n i and then after 
differentiation of ( 13 . 1 2[) with respect to n' the cross section of excitation is obtained 

da n i(vp,e) dE nl 



Vnl^n'(Vp) 

in 



T n 'T, 



ni 



de 
V2 dH(vp, 



dn' 



2[E nl - V(r) 



dr 



de 



Pni(r) 



(3.15) 



e— E n i — E n i 

where T n i is the period of oscillation over r in the final state n! . In (13. 15[) the fact was 
used that in the semiclassical approximation the derivative of the energy with respect 
to the quantum number is equal to the classical frequency: dE n i/dn' = 2tt /T' which 
follows straightforwardly from the Bohr-Sommerfeld quantization condition. 

Figure 9a shows the cross section for electron-impact ionization of helium in the 
meta-stable state He(2 3 S'). It can be seen that classical expressions give good agreement 
with experimental results in the whole range of impact energy. In Fig. 9b the comparison 
of asymptote ( I3.14p with experimental data for ionization of hydrogen atom by proton 




Figure 9. a) The ionization cross section of He( 3 S) metastable as a function of 
electron impact energy. The full curve is the classical cross section (|3.13j) and the solid 
circles are experimental data |50j . b) The ionization cross section of hydrogen atom by 
proton and He 2+ as a function of scaled impact energy. The full curve is the asymptote 
(|3.14|) . the solid and open circles are the experimental data [51] for proton and He 2+ 
as a projectile. 



and He 2+ is presented. The discrepancy at very high impact velocity in both figures 
is explained by logarithmic-like dependence of the cross section (a oc \nv p /v 2 ) due to 
the contribution from individual dipole transition. However, according to Ehrenfest's 
theorem, the classical description deals with the wavepackets. In order to obtain a 
wavepacket, which imitates the motion of a classical particle, one needs to perform 
integration over the quantum numbers in a narrow range. As a result of this averaging, 
the logarithmic energy dependence disappears. It explains the absence of the logarithmic 
energy term in the classical cross sections. 

To calculate the state selective excitation nl — > n'l', the angular- momentum transfer 
is needed. The square of the angular momentum in the final state can be written 

L' 2 = [r x (p + Ap)} 2 = L 2 + r 2 (2p ± Ap ± cos a + Ap]_), (3.16) 

where p± and Ap± are the projections of p and Ap onto the plane perpendicular to r, 
and a is the angle between p± and Ap^. Since for transition from the (nZ)-subshell the 
initial distribution q n i(r) is isotropic, the distribution with respect to L' for fixed values 
of r and b is 

, rls 1 da 2L' 
A r, b, L ) = — —— = 3.17 

dL ^J{ti-Apl){Apl-ti) 



2tt 



where = [L! ± L)/r. After averaging of (I3.17P over the impact parameter b and the 
initial distribution q n i(r) the final result is obtained as a single radial integral over the 
classically allowed region which is common to the initial and final states. The detailed 
analysis of this cross section can be found in [48l 149] . In the case n,n' 3> An and 
L, L' 3> AL (An = \n' — n\, AL = \Lf — L\) in high-energy limit the cross section has 
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the form 

32n 3 (n 2 - L 2 ) 3 / 2 

Onl^n>v{v p ) - 97rv 2( nAn 2 + nAL 2 _ 2LAnAL) 2 ( ^ 

This expression reveals some important features. For fixed values of n and L, the cross 
section sharply decreases with increasing An and AL, as (An) -4 and (AL)~ 4 . 



4. Classical representation in quantum mechanics 

At first sight, the classical representation is impossible in quantum theory, for instance, 
because of tunneling phenomenon, which is, at first sight, incompatible with the classical 
approach. However, if we consider a particle scattered by a parabolic barrier, we find 
that, on the one hand, there is an effect of sub-barrier penetration and, on the other 
hand, the time- dependent Green function, which carries full information about the 
quantum system, is exactly expressed in terms of the classical action along classically 
allowed trajectories only. Thus, in principle, the tunneling effect could be explained in 
terms of classical trajectories. 



4-1. Abel transform to the classical representation 

For the sake of transparency we consider a symmetric potential V(x) = V(—x), where 
V(0) = and V(x) increases monotonically to infinity on the semiaxis x > 0. To 
introduce the classical representation the linear equation for the quantum probability 
density p n {x) = ipn( x ) ls needed , which is obtained from the Schrodinger equation [61] 

_^f^ +2[vix) _ E /j!M + ^m Pnix)=0 . (4,i) 

This third-order differential equation has three linearly independent solutions - two of 
these are related to the Schrodinger equation and the third to the solution of the Milne 
equation having the meaning of the quantum wavelength (see Sec. 5). 

Now let us represent the quantum state n as an ensemble of classical states in the 
potential V(x) distributed in the energy e with a certain probability density 4> n (e) 



where 



is the x-distribution for the classical state at given energy e and T(e) is a classical period. 
A remarkable property of the representation H4.2[) is its reciprocity. This becomes evident 
if instead of x the potential V is taken as an independent variable. Such a change of 
variable is standard in classical mechanics and reduces the conversion of expression ( 14. 2 j) 
to the well-known Abel problem, whence 

T(e) r 1 dp n (x(V)) 



Pn(x) = q(e,x)(j) n (e)de, (4.2) 

JV(x) 



V2m 

q(e,x) = — -j== (4-3) 
T(e)Je-V(x) 



7T\ 



2m Je ^V~^~e dV 



-dV, (4.4) 
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where x(V) is the inverse of function V(x). 

The transforms (JH2J), (I4.4p are an analogue of the Fourier transform connecting 
configuration and momentum representations and it is called 'a classical representation', 
since the kernel q{e,x) is the classical probability density and 0(e) has the meaning of 
the energy distribution in the classical ensemble. 

In the classical representation, the Schrodinger equation (14. ip transforms into an 
equation for scaled energy distribution n (e) = n (e)/T(e): 

m Je a/r 



(e - E n )Me) = - I Q^ey^p-dyi, (4.5) 



with the kernel 

rX2 



. 1 ^ - V(x))V\ 
Q(/i ' £) = 15^ L T^)=e & ' ( } 



where x\ and x<i are the turning points defined by the conditions: V[x\) = e and 
V{x2) = yU. The discrete spectrum E n is obtained from the boundary condition: 
n (e) — > as e — > oo. Equation (14. 5 p looks like a balance equation for 'virtual' transitions 
between the classical states with the energies e and fi. It is very interesting to extract 
the probability of such transitions from equation (14. 5p . 

The energy distribution n (e) obeys the obvious normalization condition 







<Pn{t)dt = 1 (4.7) 

In addition to this, the general condition holds 

roo 

/ e0 n (e)de = E n , (4.8) 
Jo 

which is not trivial and means that the average energy of classical ensemble is equal to 
the eigen-energy of quantum state \n). 

4-2. Harmonic oscillator 

Let us consider the harmonic-oscillator problem 

V(x) = \nu?x 2 (4.9) 

as an example in which the classical representation is employed. In this case the kernel 
can be calculated explicitly: 

Q(^e) = ^moo 2 (fi-2e). (4.10) 

Substituting this into (14. 5p yields an equation for harmonic oscillator in classical 
representation 

ftV r d 2 <j) n {e) dcf) n (e) 
4 L de 2 de 

which formally coincides with the radial Schrodinger equation for the Sturmian basis 
functions of the hydrogen atom with the negative 'angular quantum number' / = —1/2. 



+ (E n - e)0 n (e) = 0, (4.11) 
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In the Sturmian approach the energy is fixed and the charge is quantized. The energy 
E n plays the role of the 'charge' in this 'hydrogen Schrodinger equation', its spectrum 
is E n = hu]{n + 1/2), and the corresponding energy distribution reads 

= ^e- 2 ^L n (4e/M, (4-12) 
where L n (x) is a Laguerre polynomial. 

4-3. Inelastic transitions in the Feynman model 

As practical application of the classical representation, the problem of description of the 
initial and final states in classical trajectory Monte Carlo (CTMC) method [52] in the 
theory of atomic collisions is the first candidate. In this method the initial state of the 
electron is specified as an ensemble of classical states with the energy equal to the energy 
E n of the corresponding atomic state. Then the evolution of this ensemble is calculated 
according to the classical equations of motion. The main difficulty in this approach is to 
extract the contribution of different quantum states in the final ensemble as t — > oo, since 
the final energies of classical trajectories assume continuous values that do not coincide 
with the atomic spectrum. On the other hand, the classical representation (14 .4p provides 
the exact description of quantum states in terms of an ensemble of classical trajectories. 
This fact can be used to modify the CTMC method in the treatment of the initial and 
final states. Of course, the solution of the dynamical problem remains classical, i.e. 
approximate. To illustrate such an approach, let us consider inelastic transitions in the 
Feynman model (63] with the nonstationary potential 

V{x, t) = \nu?x 2 + a(t)x, (4.13) 

where a(t) is the strength of a homogeneous external field that is time-dependent and 
tends to zero as t — > ±oo. In this case, the general solution of the classical equation of 
motion can be obtained explicitly 



I /"* v2e 
x(t) = / a(t') sin cu( t - t')dt' + - — smuCt-r), (4.14) 

mu J-oo uj 

where e and r are the initial energy and phase. As t — > +oo this solution describes 
harmonic oscillation with an energy 

/i = e + v + 2\fev cos r, (4. 15) 

where 

1 r+oa 2 

Since the distribution over the initial phase is uniform, the energy distribution in the 
final state at initial energy e has the form 

fa, e) = ~ = y-^fo-/*) , (4.17) 
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where /ii = (y/e — \/v) 2 , /i 2 = \/^) 2 and 0(x) is the unite step-function. Then the 

probability of an inelastic transition from the state \n) to the state \k) is [61] 

P nk (l) = [°° drfM f 2 dep{p, e)cj> n {e) = ^(7)^(7))^, (4.18) 

where <pk{p) and 0„(e) are the energy distributions of the initial and final states 
defined by fl4.12j) . e x = - a/z^) 2 , e 2 = (y/Jl + v^) 2 > 7 = vjfujj and £$(x) is 
a generalized Laguerre polynomial. The probability (14. 18j) coincides with the exact 
quantum expression (see [63J ) , although the classical solution (I4.14p is approximate 
from a quantum point of view. 



4-4- Semiclassical approach 



The rough classical approximation is obtained if we neglect the right-hand side in 
equation (14.51) . Then <p n (e) = S(e — E n ) and from (14.21) the density of probability 
takes the form 

. , J 2m 

p n (x) = (4.19) 



(4.20) 



T{E n ) x jE n -V{x) 

However, the solution (j) n (e) as well ) has essential singularity as h — > 0: 



^ocsin^^e)). 



To construct semiclassical approximation, the integral in the right-hand side of equation 
( 14. 5 p has to be integrated by parts twice 



h 2 


rdV(x) 




2d 2 (t) n {e) 
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dx 


V=e- 


de 2 



+ 



h 2 d 2 v(x) 



m dx 2 



v= 



de m 



d 2 Q(fi,e) d4> n (/j,) 



dp 2 



d/j, 



dfi. 



(4.21) 



Here the coefficients in front of the second and first derivatives of n (e) are obtained 
from definition of Q(fj,,e) (14. 6p . Since </>„(e) has essential singularity with respect to h, 
in zero order only the first term in the right-hand side should be retained. Then W(e) 
takes the form 



W(e) 



2m[E n 



vdV(x) 




dx 


V=e>. 



-1 



de'. 



(4.22) 



After substitution P^ZUj) with W(e) (|4T2^1) into expression (T4T2]) . the density of 
probability can be estimated by the expression 



p(x) oc f g(e,x)sin(| / J2m(E n - V) 
Jv(x) K a J v 



dV{x) 



dx 



V=e' 



de')de. (4.23) 



Semiclassical expansion of p(x) is generated by consecutive integration by parts in (I4.23p . 
In the leading order 

-2 



(4.24) 



p(x) oc cos(- J X ^2m(E n - V{x'))dx') = 
= 2 cos 2 (- J X \j2m{E n - V(x'))dx'/h) - 1 . 
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The last term, -1, has to be omitted because it is beyond the semiclassical series, and 
the result coincides with the well-known semiclassical expansion of p(x) in the leading 
order. 

4-5. Scattering problem; tunneling phenomenon 

In quantum mechanics, the scattering problem is out of Hilbert space and we have no 
well-defined physical stationary states because of normalization problem. It is a very 
serious defect of the theory. Below the boundary of the continuum, the aim of the theory 
is the search for stable (or stationary) states which exist at discrete energy values only. 
Above the boundary of the continuum the classical motion of particle is, in principle, 
not stable (the particle comes from infinity and then goes to infinity as t — > oo), and 
we cannot clearly formulate a physical problem like in the case of bound states. But, 
nevertheless, if we accept the main idea of classical representation - the quantum state 
is reproduced by an ensemble of classical states - then we can interpret, for example, the 
tunneling phenomenon as penetration through a barrier of part of this ensemble with 
energy of trajectories above the top of a barrier. 

5. Semiclassical series 

5.1. Higher orders of semiclassical series for quantum wavelength 

The semiclassical approximation for the one-dimensional Schrodinger equation 



is the asymptotic expansion of the wave function ip(x) with respect to the small Plank 
constant H. Within the standard semiclassical approach, when the logarithm of the wave 
function is expanded 



the treatment of higher orders encounters difficulties because the recurrence relation, 
from which the n-th term S n is obtained, involves in a nonlinear manner all previous 
terms d\, k < n [HZ]. To avoid this difficulty, the Milne transformation should be used, 
which reduces the problem to the study of semiclassical expansion of the 'wavelength' 
A(x), related to the wave function by 




(5.1) 



oo 



hii/>(x) = Yl 



(5.2) 



n=-l 




(5.3) 



and satisfies the third order linear differential equation (see also Sec. 4.1) [61] 




(5.4) 



In this approach the coefficients of the expansion 



oo 



A(x) = ][> 2n A n (x) 



(5.5) 



n=0 
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are subject to the two-term linear recurrence relation 

w v 1 f x dx' d 3 A n _i(x') C n 

X n (x) = / — h , (5.6) 

/p(i) J •Jpix') ax yp( x ) 



where = y2m[E — V{x)\ is classical momentum and \q{x) = l/p(x). The constant 
C n is determined by a nonlinear Milne equation; the whole effect of constants C n in 
recurrence relation (15.61) . as well as in solution of equation (I5.10p . is that their change 
results only in the change of the normalization factor (in detail see [61]) and, without 
loss of generality, they will further be omitted. 

Now let us consider the divergence of the semiclassical expansion for a particle with 
zero energy in a power-law potential taken in the form V(x) = — a 2 x 2u /2m, where a 
and v are arbitrary constants. This is an interesting case in that it describes all types 
of violations of the semiclassical approximation because of power-law singularities. The 
value v — 1/2 corresponds to an ordinary turning point. From the recurrence relation 
(I5.6P it is easy to find the general term 

-i) n r(n + ^g)I> + |)r(n + ^±%) , „ + 1 



K{X) ~ ^r(^)r(i)r(^) U, + J *-"> ( 5 - 7 ) 

From this expression follows that for v < — 1 the semiclassical expansion breaks down 
as x — > oo, while for v > —1 it breaks down as x — > 0; correspondingly, it becomes 
exact in the limit x — > in the first case and in the limit x — > oo in the second. When 
v — — 1± l/(2g + 1) (q = 0, 1, 2, ...), the expansion is cut off at n = q, and summation of 
the first q terms gives the exact solution of the Schrodinger equation in the form (15. 3j) . 

Expansion (15. 7p diverges as (n!) 2 , so that the corrections A n can be taken into 
account only until they begin to grow. The index N of the order at which the expansion 
should be truncated is determined by setting the derivative of A n with respect to n equal 
to zero. At large n the following result is obtained N = ax u+1 /K(y + 1), which can be 
rewritten as 

N — — p{x')dx', (5.8) 

where xq is the point where the semiclassical approximation breaks down (xq = if 
v > —1 or xq = oo if v < —1). Thus, the critical index N is equal (in terms of h) to 
the classical action from the point x to the singularity xq] it is invariant in form and is 
found for arbitrary power-law singularities. 

5.2. Renormgroup symmetry 

The semiclassical expansion can be written, as well, in the form [651 ES] 

^(x) = 1 e ^(»)/»f;[±£ ] n Mx ), (5.9) 

where S(x) is the classical action calculated from the singular point xq. Further, it 
is convenient to take the action S(x) as a new independent variable instead of the 
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coordinate x. Then, the substitution of (15. 9p into the Schrodinger equation leads to the 
recurrence relation 

where 

^-jM- (511) 

and the leading term is <po(x) = 1. Since Q(x) has a pole of second order at the singular 
point xq with residual v{y + 2)/[2(i/+ l)] 2 , the solution of the recurrence relation (15.101) 
has a pole of n-th order with factorial growing residual: 

MS) « (-i)"^^- ( 5 - 12 ) 

To analyze higher-order terms, let us introduce asymptotic expansion of <p n (S) in inverse 
powers of n in the following special form: 

«S) = (-l)"[^tf'(S) + I^tf>(S)... 



- + £ ^4 1) (s) + ...l 



with the initial condition (f)^\x) = 1. Then, from equation (I5.10p . for (^^(S) the 
recurrence relation is obtained 

d$\S) d^Us) (1) 

— ^ — = — Jcp — + Q( s ) ( t ) k-i(S), (5.14) 

which is exactly the same as (15.101) . Moreover, if one expands the function (p^\S) for 
large k in the same manner 



- + £ §^ £l #^) + 4 ( 5 - 15 ) 

then the new function <p^\S) will again be subject to the recurrence relation (I5.10p . 
We can repeat this procedure m times and obtain the same recurrence relation, i.e., the 
same sequence {(j)^ 1 ' (S)} as {(p n (S)}. This amazing property reminds a self-similarity 
structure appearing in the theory of dynamical systems where the picture of irregular 
(chaotic) motion repeats itself as one goes to finer scales of phase space - so called 
renormgroup symmetry [50] ■ This property was discovered by Dingle (see [B3], Chapter 
XIII) and independently in [66] . 



5. 3. Criterion of accuracy 

With increasing n the terms of the asymptotic expansion (15.91) at first decrease until for 
some number N they reach the minimal value, after which they increase. The accuracy 
of the finite sum is determined by the last retained term. Therefore, it is correct to sum 
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up only first N terms; the addition of the higher terms would make the result worse. As 
in Sec. 5. 4.1, the critical index N is determined by setting the derivative of h n <ft n (x) with 
respect to n equal to zero. At large values of n, using (15. 12ft and Stirling approximation 
for the n\, it follows 

iV=-r-. (5.16) 
n 

i.e., the classical action up to nearest (in terms of S) singularity determines the number 
of terms of expansion (15.9p for best accuracy. Taking into account that the expansion 
in Sec. 5. 4.1 is over even powers of h, one can see that (15.81) and (15.161) coincide, but the 
result (15.161) is obtained in the general case. The relative error is estimated by 

(tr^TT ~ < 5 - i7 » 

It is also interesting to see how many terms it is allowed to keep in asymptotic sub- 
expansion (15.13!) . For large k we can replace (x) by the leading term of expansion 
(I5.15p . Then, it can be verified in the same manner as for derivation (I5.16P that the 
minimal term in sub-expansion (15. 131) is at k = n/2. The same result is obtained for all 
next sub-expansions (j)^(x) (m = 2,3,4, ...). 



6. Concluding remarks; classical physics and measurement 

We obtain any data from measurement which, in principle, gives the result in terms of 
classical physics. Intuitively, it is clear what the concept 'measurement' means; however, 
it cannot be formalized, i.e. to be written in mathematical form. What is the relation 
between quantum theory and measurement? In classical mechanics and classical electro- 
dynamics the concepts 'material points' and 'electro-magnetic waves' were introduced. 
In this sense we can consider quantum theory as a theory of 'information field' which is 
described by wave function \I/(r) [68J; it is additional to 'material points' and 'electro- 
magnetic waves' substances but without mass and energy. The problem of the existence 
of 'information field' is similar to the situation with electro-magnetic field - nobody 
suspected its existence before classical electro-dynamics. Experimentum crucis, which 
proves the existence of information field, is based on the Einstein-Rozen-Podolsky 
phenomenon. One of the realizations of this experiment is the following. After a 
radiative decay of a hydrogen atom from a metastable 25'-state two photons are emitted 
having spins with opposite directions 

s 1 = -s 2 , (6.1) 

since the initial (25*) and final (IS*) states of hydrogen atom have angular momentum 
equal to zero. However, the orientation of each spin is uncertain in the same manner 
as a position of a particle in a potential well. If the measurement of the direction 
of the spin of the first photon is performed the second photon takes the opposite 
direction of spin, according to (16. ip . independently of the distance between photons. 
The actual fixation of the spin orientation of the second photon happens immediately 
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and it is an actual changing of its state (Bell's inequality [HH]). It does not contradict 
the relativistic restriction v < c because the carrier of information is not material, i.e. it 
has no relation to the transfer of mass or energy. Recent experiment [70] demonstrates 
that the speed of quantum information is at least 10 4 greater than the speed of light 
c. This experiment confirms that an information field does exist. Notice that this 
interpretation is in contradiction with wave-particle duality: the wave-type Schrodinger 
equation describes just 'information field' ^(r). Thus, there is no duality because of two 
different physical substances - particles and information field. The interaction between 
particles and information fields is measurement which cannot be formalized, i.e. to be 
written in mathematical form. 
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